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Voltage-gated potassium channels elicit membrane hyperpolarization through voltage-sensor 
domains that regulate the conductive status of the pore domain. To better understand the 
inherent basis for the open-closed equilibrium in these channels, we undertook an atomistic 
scan using synthetic fluorinated derivatives of aromatic residues previously implicated in the 
gating of Shaker potassium channels. Here we show that stepwise dispersion of the negative 
electrostatic surface potential of only one site, Phe481, stabilizes the channel open state. 
Furthermore, these data suggest that this apparent stabilization is the consequence of the 
amelioration of an inherently repulsive open-state interaction between the partial negative 
charge on the face of Phe481 and a highly co-evolved acidic side chain, Glu395, and this 
interaction is potentially modulated through the Tyr485 hydroxyl. We propose that the 
intrinsic open-state destabilization via aromatic repulsion represents a new mechanism by 
which ion channels, and likely other proteins, fine-tune conformational equilibria. 
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Voltage -gated potassium (Kv) channels are membrane 
proteins whose activity underlies the regulation of basal 
electrical excitability and action potential repolarization in 
virtually all electrically excitable cells 1 . Significant progress over 
the past decade in the structural biology of membrane proteins 
has advanced our understanding of ion channels (and especially 
Kv channels), likely to a much greater extent than any other 
membrane protein family 2 ' 3 . It is now well established that these 
proteins are formed by four subunits, arranged symmetrically 
around a central, ion-conducting pore. Each subunit contains six 
transmembrane domains (S1-S6), with S1-S4 comprising the so- 
called voltage-sensor domain and S5 and S6 making up the pore 
domain. Under physiological conditions, Kv channels cycle 
rapidly between open and closed states and the relative stability 
of these states is of critical importance for a number of 
physiological processes 4-6 . 

The conductive state of the pore represents the functional 
summation of multiple input energies, including the inherent 
equilibrium between open and closed conformations of the pore, 
as well as the 'work' that is performed by the voltage- sensing 
domains to bias this equilibrium. In terms of the former, 
computational and functional approaches support a model 
whereby the closed, nonconductive conformation of the pore is 
more stable compared with the open conformation, or conversely, 
the open pore represents a meta- stable complex 7-9 . The basis for 
this relative open-state instability is poorly understood, but 
numerous approaches have converged on the importance of side 
chains near the S4-S5 linker and the distal S6 segment, which 
line the interface between the voltage-sensing domain and 
cytoplasmic face of the pore region 10 . Interestingly, many of 
these side chains are aromatic in nature and are often highly 
intolerant to conventional site- directed mutagenesis. 

We turned to the in vivo nonsense suppression method for 
incorporation of unnatural amino-acid side chains 22 , to subtly 
alter the aromatic side chains in an atom-by-atom manner, thus 
allowing for a more careful dissection of their physico-chemical 
contributions. The results show that reducing the negative surface 
potential on the face of a single highly conserved aromatic side 
chain, Phe481, through serial fluorination, leads to an increase in 
open-state stability, as apparent from conductance-voltage (GV) 
relationships, channel closing kinetics and kinetic modelling. By 
contrast, fluorine substitution at four nearby aromatic side chains 
(Phe401, Phe401, Phe484 and Tyr485) has little functional 
consequence. Interestingly, neutralization of a highly co-evolved 
nearby acidic side chain, Glu395, abolishes the effects of 
fluorination of Phe481, an effect not observed when 
neighbouring basic side chains are neutralized. Taken together, 
these results demonstrate for the first time the functional 
importance of the electrostatic surface potential of an aromatic 
side chain for biasing the conformational equilibria of an ion 
channel. Further, they are consistent with the presence of an 
intrinsically repulsive electrostatic interaction between the face of 
an aromatic and the net negative charge of a highly co- evolved 
acidic side chain (located on an adjacent subunit) that may serve 
to inherently destabilize the channel open state, an interaction 
potentially modulated through the Tyr485 hydroxyl moiety. 

Results 

Phe481 electrostatic surface potential regulates function. 

Structural and functional data have highlighted the importance of 
five aromatic side chains in the cytoplasmic end of S5 and the 
distal part of S6 (Fig. la). Importantly, four of these aromatic side 
chains (Phe401 and Phe402 in S5, Phe484 and Tyr485 in S6), 
have been suggested to have an important role in channel func- 
tion 10-13,19 . Further, mutations of Phe481 in S6 to nonaromatic 



side chains such as Ala, Asp , Leu or Cys (Supplementary Fig. 
SI) result in nonfunctional channels and even the most 
conservative conventional mutation (Phe481Tyr) results in 
significant functional perturbations (Supplementary Fig. SI). 
Here, we sought to more clearly define the role of these five 
aromatic side chains in channel function, as interpretation of 
results obtained by conventional site-directed mutagenesis is 
hampered by the fact that mutations can change many physico- 
chemical properties (such as size, hydrophobicity, H-bonding 
ability or charge) of a side chain. To more carefully characterize 
the contribution of these five aromatic side chains to pore 
stability, we turned to the in vivo nonsense suppression method 
for incorporation of unnatural amino acids 22 . In particular, we set 
out to assess the contribution of the negative electrostatic surface 
potential of these aromatic side chains to the inherent open 
versus closed equilibrium by using fluorinated Phe derivatives. 
These valuable tools specifically reduce the negative electrostatic 
surface potential of the Phe side chain, while leaving size, shape 
and hydrophobicity virtually unchanged. Fluorinated Phe 
derivatives have been extensively used to identify cation-pi 
interactions, a predominantly electrostatic interaction between 
the electronegative surface potential of an aromatic side chain and 
an organic cation 24-26 , as the addition of each fluorine leads to a 
stepwise decrease in the electronegative surface potential. To test 
if channel function is sensitive to changes in the electrostatic 
surface potential of Phe401, Phe402, Phe481, Phe484 or Tyr485, 
we incorporated F 3 -Phe at each position, a manipulation that 
effectively abolishes the negative surface potential on the aromatic 
face. Surprisingly, trifluorination had no measurable effect on 
the GVs at positions 402, 484 and 485 and only produced a 
modest (<10mV) right-shift in position 401 (Fig. lb-e). In 
contrast, introducing F 3 -Phe at position 481 led to a significant 
(P<0.01) hyperpolarizing shift (~30mV) in the GV (Fig. If, 
Supplementary Table SI), a particularly noteworthy phenotype 
because of its location at the critical interface of S6 and the S4-S5 
linker 14-21 ' 27 ' 28 and its high conservation among Kv channels 
(Fig. 2a,b). Importantly, alterations of channel gating produced by 
F 3 -Phe substitution in Phe481TAG channels was not due to 
nonspecific effects evoked by the in vivo nonsense suppression 
method as re-introducing the native Phe by nonsense suppression 
produced channels that were functionally indistinguishable from 
their wild-type (WT) counterparts (Supplementary Fig. S2, 
Supplementary Table SI). In contrast, replacing Phe481 by 
cyclohexylalanine, a cyclic but nonaromatic derivative of Phe 
likely resulted in nonfunctional channels, as cyclohexylalanine- 
injected oocytes failed to display ionic currents (Supplementary 
Fig. S2). These results demonstrate that a cyclic side chain alone is 
not sufficient for normal channel function and strengthen the 
hypothesis that the electrostatic surface potential of Phe 481 has a 
critical role in channel function. 



Phe481 mediates open-state stability. Although the incorpora- 
tion of F 3 -Phe at position 481 resulted in a marked functional 
phenotype, confirmation of a contribution of the negative elec- 
trostatic surface potential of the aromatic face requires demon- 
stration of additive energetic effects of serial fluorination of the 
Phe side chain. This is accomplished by introduction of mono- 
and difluorinated Phe derivatives (F-Phe, F 2 -Phe), which dimin- 
ish the negative electrostatic surface potential of Phe481 in a 
stepwise manner. As shown in Fig. 2c, each added fluorine pro- 
duced a progressive left-shift in the GV and a steeper voltage- 
dependence (Supplementary Table SI) and the free energy 
differences between open and closed states for WT and mutant 
channels (ZFV 1/2 7 ) correlated well with the number of added 
fluorines (Fig. 2d, Supplementary Table SI). The activation rates 
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Figure 1 | The effect of fluorination is specific to position 481. (a) Model of the pore domain of the tetrameric Kv1. 2/2.1 chimera (the voltage-sensing 
domains are omitted for clarity), with the positions of Phe401, Phe402, Phe481, Phe484 and Tyr485 highlighted in yellow. The inset shows a close-up 
view of the five aromatic side chains; (b-f) GVs for WTand F 3 -Phe incorporated at positions 484 (b), 485 (c), 401 (d), 402 (e) and 481 (f); note that 
at position 485 the native Tyr was also replaced by Phe; insets show currents ( - 60 to +30mV in (a,c-f); -60 to +90mV in (b); all in 10 mV 
increments) and electrostatic surface potential maps for Phe and F 3 -Phe (ESP, red = — 15 kcal mol _1 , green = 0 kcal mol _1 , blue= + 15 kcal mol _1 ). 
Scale bar (horizontal), 20 ms; Scale bar (vertical), 5|iA in (c,d-f), 2|iA in (c). WT, V 1/2 = - 23.0 ± 0.5 mV; Phe401TAG + F 3 -Phe, V 1/2 = - 14.4 ±0.7; 
Phe402TAG + F 3 -Phe, V 1/2 = - 19.0 ± 0.6 mV, Phe484TAG + F 3 -Phe, V 1/2 = - 24.6 ± 0.6 mV; Tyr485TAG + Phe, V 1/2 = 1.9 ± 5.1 mV, Tyr485TAG + 
F 3 -Phe, V 1/2 = 4.6±3.9mV; Phe481TAG + F 3 -Phe, V 1/2 = - 45.7 ± 0.8 mV (n = 4-10 each). All data = mean ± s.e.m. 



of the ionic conductance were virtually identical for all constructs number of substituents (Supplementary Fig. S3). The 
(Fig. 2e), whereas the deactivation rates were slowed significantly fluorination-induced left-shift in the GV, concomitant with the 
with each added fluorine (Fig. 2f) and correlated well with the significant slowing of the deactivation rates (and unchanged 
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Figure 2 | Reducing the negative electrostatic surface potential on Phe481 stabilizes the channel open state, (a) Model showing the tetrameric 
Kvl. 2/2.1 chimera from the intracellular side. The distal part of S6 has been removed for better visibility of Phe481 (Shaker numbering), which is highlighted 
in yellow, (b) Sequence alignment of the distal end of S6 of various Kv channels: Shaker (Gl:288442), Kvl.l (Gl:119395748), Kv2.1 (Gl:84570020), 
Kv3.1 (Gl:298603), Kv4.1 (Gl:8272404), Kv5.1 (Gl:24418476), Kv6.1 (Gl:24418479), Kv7.1 (Gl:6166005) and Kv8.1 (GI7657289). The conserved 
Phe (Phe481 in Shaker) is highlighted in yellow, (c) GVs of Phe and fluorinated Phe derivatives incorporated at position 481 (n = 8-10 each; see also 
Supplementary Table S1); insets show representative currents (-70 to +40mV in 10 mV increments, F 3 -Phe reproduced from Fig. If) and electrostatic 
surface potential maps for Phe and F 3 -Phe (scale as in Fig. 1). (d) ZFV 1/2 values plotted against the number of fluorines for Phe and fluorinated Phe 
derivatives incorporated at position 481. (e) Activation rates plotted against voltage for Phe and fluorinated Phe derivatives incorporated at position 481; 
insets show sample traces for Phe and F 3 -Phe (depolarizations in 10 mV increments from a holding potential of -80 mV, traces shown from — 70 mV to 
40 mV). (f) Deactivation rates plotted against voltage for Phe and fluorinated Phe derivatives incorporated at position 481; insets show sample traces 
for Phe and F 3 -Phe (after a depolarizing step to 0 mV different test potentials were applied in 10 mV increments, traces shown from 0 to -100mV). 
Scale bar (horizontal), 20 ms in (c,e), Scale bar, 50 ms in (f); Scale bar (vertical), 10 uA in (c,e), Scale bar, 5u.A in (f). All data = mean ± s.e.m. 



activation rates) suggested that stepwise fluorination progressively 
stabilizes the channel open state. 

To test the validity of this hypothesis, we turned to kinetic 
modelling based on the Zagotta-Hoshi-Aldrich model of Shaker 
channel activation 29 (Supplementary Fig. S4). We systematically 
simulated the effects of altering individual rates in the activation 
pathway. The simplest manipulation that recapitulates our 
experimental findings is to simply decelerate the closing rate for 
the final 'concerted' gating step (highlighted by a green box in 
Supplementary Fig. S4). This alteration of the model generates a 
leftward shift and steepening of the GV relationship, and slower 
deactivation rates that closely match the experimental effects of 
fluorination of Phe481 (Fig. 3a,b, see Supplementary Fig. S4 for 
specific rates and valences of model transitions). Perturbation 
of other model parameters failed to generate this combination 
of effects. 

If, indeed, fluorination affected only the final opening 
transition, then the fluorination-induced phenotype should be 
retained in any mutation that separates the final opening 
transition from any of the earlier transitions, such as those 
identified by Aldrich and co-workers 30 ' 31 in the S4 segment of 
Shaker channels (Val369Ile, Ile372Leu, Ser376Thr or ILT). We 
introduced a subset of these mutations (Ile372Leu, Ser376Thr or 
LT) on the Phe481TAG background and introduced either Phe or 
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F 3 -Phe at position 481. As shown in Fig. 3c, simultaneous 
introduction of the Ile372Leu and the Ser376Thr mutation (LT) 
resulted in a >80mV right-shift in the GV. However, 
trifluorination left-shifted the GV by 55 mV towards WT and, 
importantly, the free energy difference between Phe and F 3 -Phe at 
position 481 was virtually identical on both the WT and the 
mutant (LT) background (Fig. 3d and Supplementary Table S2), 
strongly suggesting that fluorination does, indeed, stabilize the 
channel open state by altering the open-to-closed rate of the 
concerted pore gating step. 

Effects of fluorination on S4 movement and inactivation. 

Despite its location in the distal part of S6, Phe481 also resides in 
close proximity to the S4-S5 linker and could thus directly affect 
electromechanical coupling and/or voltage-sensor movement. 
Further, although C-type (or slow) inactivation is mediated by 
side chains near the selectivity filter 32 , mutations distant from this 
domain have been shown to affect the process of C-type 
inactivation. Consequently, we set out to test directly if 
fluorination of Phe481 has effects on S4 movement and/or 
C-type inactivation. To address the former, we used voltage- 
clamp fluorometry (VCF) 33 ' 34 to track S4 movement with a 
fluorophore as a reporter. We found that labelling Ala359Cys 
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Figure 3 | Kinetic modeling and mutational studies confirm open-state stabilization through fluorination. (a) Experimentally observed GV data points 
for Phe and F 3 -Phe (open and blue circles, respectively; reproduced from Fig. If) overlaid by the simulated data based on kinetic modelling using the ZHA 
model (red lines); (b) Experimentally observed deactivation rates for Phe and F 3 -Phe (open and blue circles, respectively; reproduced from Fig. 2f) 
compared with the simulated deactivation rates for Phe and F 3 -Phe based on kinetic modelling using the ZHA model (red lines); (c) GVs for Phe and F 3 -Phe 
introduced on the lle372Leu, Ser376Thr (LT) double-mutant background (open and blue circles, respectively). The fit for the WT GV is shown in grey, 
insets show representative currents. Scale bar (horizontal), 20 ms; Scale bar (vertical), 4 uA (d) Bar diagram comparing the free energy difference between 
Phe and F 3 -Phe incorporated at position 481 on the WT (empty bar) or LT double-mutant background (black bar). All experimental data = mean ± s.e.m. 



(located at the top of S4) with tetramethylrhodamine-5- 
maleimide with either Phe or F 3 -Phe in position 481 resulted 
in robust fluorescence signals that displayed an identical 
voltage- dependence (Supplementary Fig. S5), suggesting that 
fluorination of Phe 481 specifically alters the inherent open versus 
closed equilibrium of the pore domain but not the voltage- 
dependence of S4 movement. Next, we set out to investigate the 
possibility of functional coupling between Phe481 and the 
selectivity filter region. However, comparing the time course 
and extent of C-type inactivation with either Phe or F 3 -Phe in 
position 481 yielded identical results (Supplementary Fig. S6), 
demonstrating that fluorination does not affect C-type 
inactivation. 



Identifying potential interaction partners of Phe481. The 

functional impact of fluorination is highly specific for position 
481 and originates from changes in the electrostatic surface 
potential. To identify possible interaction partners that may be 
involved in this electrostatic effect, we turned to statistical 
coupling analysis (SCA) 35 ' 36 , a method for identification of 
co-evolved pairs of side chains (or domains) 37 (Supplementary 
Fig. S7). In particular, we focused on charged side chains within 
13 A of Phe481: four basic side chains (Lys380, Arg387, Lys390 
and Arg394) and a single acidic side chain (Glu395). In effect, we 
sought to determine whether Phe481 engages in an energetically 
attractive (cation -pi) interaction with a nearby basic side chain or 
in an energetically repulsive interaction mediated by an acidic 



side chain. Figure 4a shows a pseudo-heat map of SCA scores for 
Phe481 and all five charged side chains, illustrating the relatively 
low values for the basic residues, and surprisingly high co- 
evolution scores for Phe481 and Glu395 when compared against 
the entire channel sequence. We generated a histogram 
illustrating the frequency of SCA scores among all pairwise 
comparisons in the channel, highlighting the SCA scores of all 
five charged side chains (Lys380, Arg387, Lys390, Arg394 and 
Glu395), with Phe481 (Fig. 4b). This analysis yields a remarkable 
disparity between the four basic side chains and the single acidic 
side chain: while Arg387, Lys380, Arg394 and Lys390 all had low 
scores (0.13 kT*, 0.20 kT*, 0.21 kT* and 0.47 kT*, respectively, in 
terms of 'statistical energy'), Glu395 had one of the highest scores 
in the entire data set (1.93 kT*) (Fig. 4b), indicating a particularly 
high degree of co-evolution between Phe481 and Glu395, an 
observation further highlighted by the close physical proximity of 
the two side chains (see structural model in Fig. 4b). Additional 
evidence for a potential interaction between Phe481 and Glu395 
stems from the observation that in Kv sequences with a Phe in the 
position equivalent to Shaker 481, only the position equivalent to 
Glu395 shows a strong preference for a single, acidic side chain 
(Glu), while positions 380, 387, 390 and 394 show little or no side 
chain preference (Supplementary Fig. S8). 



Phe481 does not interact with nearby basic side chains. Our 

functional data show that Phe481 is likely engaged in an elec- 
trostatic interaction, and SCA suggests that four nearby basic side 
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Figure 4 | SCA analysis highlights the varying levels of co-evolution between Phe481 and nearby charged side chains, (a) SCA scores of Phe481 and all 
five charged side chains that are within 12 A of Phe481 with the vertical axis representing all side chains considered for the SCA analysis (corresponding 
to positions 208 to 487 in Shaker, with black boxes representing a-helical regions in the channel). Scale bar as in (b). (b) Bar diagram indicating the 
percentage of SCA scores exceeding the binned threshold indicated on the x axis for all pairings in the SCA analysis. The SCA scores of all five charged 
side chains with respect to Phe481 are highlighted with arrows (Arg387 = 0.13 kT*, Lys380 = 0.20 kT*, Arg394 = 0.21 kT*, Lys390 = 0.47 kT* and 
Glu395 = 1.93 kT*); Note the high co-evolution score for the Glu395-Phe481 coupling. Inset shows the scale of the heat map used in (a,b) and a model 
displaying all five charged side chains that are within 13 A of Phe481. Note that the charged side chains are coloured by their SCA score values (Phe481 is 
shown in dark grey as all SCA scores are relative to Phe481 in this analysis). 



chains are improbable candidates for a (co-evolved) interaction 
with Phe481. To validate our statistical approach we reasoned 
that neutralizing any of the four basic side chains should have no 
impact on the fluorination-induced phenotype at Phe481. Indeed, 
we found that neutralization of Lys380, Arg387, Lys390 or 
Arg394 had no effect on the fluorination of Phe481: we still 
observed a stepwise left-shift in the GV with fluorination 
(Fig. 5a-d) and plotting the respective AZFV 1/2 values against the 
number of fluorines resulted in a linear trend, comparable to that 
obtained with the WT background (Fig. 5f, Supplementary Table 
SI). Additionally, we investigated a possible role of Arg3 77 
because despite its large (>18A) distance to Phe481, Arg3 77, 
similar to Phe481, is highly sensitive to modest conventional 
mutations 38 ' 39 . The charge- neutralizing mutation Arg377Gln is 
nonfunctional 38 ' 39 , but if a specific electrostatic interaction 
between Arg3 77 and Phe481 was present, even replacing Arg 
by Lys should result in altered interaction energetics 40 ' 4 . 
However, the resulting free energy difference on the Arg377Lys 
background was identical to that of WT, ruling out an 
electrostatic interaction between Arg377 and Phe481 (Fig. 5e,f). 

Neutralizing Glu 395 abolishes the effect of fluorination. If the 

electrostatic surface potential of Phe481 is crucial for open-state 
stability, but nearby positively charged side chains do not appear 
to have a role, what other side chains could be involved in the 
electrostatic interaction with Phe481? Interestingly, Glu395 is not 
only highly conserved among Kv channels (Fig. 6a), it is also the 
only other charged side chain in the vicinity of Phe481 (around 
7 A away in the Kvl.2/Kv2.1 chimera open channel crystal 
structure 3 and around 5 A away in a refined model of the Kvl.2 
crystal structure 42 ) (Fig. 6b). More importantly, our SCA analysis 
identified Glu395 as the candidate with the highest co- evolution 
score of all charged residues within 13 A of Phe481 (see Fig. 4). 
Figure 6c shows the lack of a stepwise left- shift in the GV through 



fluorination at position 481 on the Glu395Gln background, 
exemplified by the incorporation of F 3 -Phe, which results in a GV 
that is indistinguishable from that of Phe incorporated in position 
481 on the Glu395Gln background (Supplementary Table SI). 
The effect is made clear by plotting the free energies of these three 
mutants, as well as mutants obtained through fluorination on a 
different hydrophilic background (Glu395Asn, Fig. 6d,e) and 
various hydrophobic backgrounds (Glu395Ala, Glu395Leu and 
Glu395Ile, Supplementary Fig. S9 and Supplementary Table SI), 
none of which show a trend with fluorination. 

A possible role for Tyr485 in modulating the interaction. Given 
that some structural data 3 suggest the Tyr hydroxyl of Tyr485 to 
reside between Phe481 and Glu395, we tested the effect of 
fluorination in position 481 on the Tyr485Phe background to 
examine a potential modulating role of the Tyr hydroxyl in the 
repulsive interaction between Phe481 and Glu395. Removal of the 
Tyr485 hydroxyl resulted in a large right- shift of the GV, as well 
as a rapid and irreversible inactivation phenotype (Supplementary 
Table SI and Supplementary Fig. S10); and fluorination in 
position 481 did have no further effect on channel function, 
potentially suggesting that the Tyr485 hydroxyl modulates an 
indirect repulsive interaction between Phe481 and Glu395. 
However, the drastic phenotype of Tyr485Phe alone could be 
indicative of major structural perturbations that may result in a 
nonspecific disruption of the interaction between Phe481 and 
Glu395, thus complicating the interpretation of the findings with 
Tyr485Phe. 

Further insight into a potential repulsive interaction. Given 
that our findings point towards an important role of the negative 
electrostatic surface potential of Phe48 1 and statistical, structural 
and functional data support the notion of Glu395 being a 
potential (direct or indirect) interaction partner, we turned to ab 
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Figure 5 | Nearby positive charges are unlikely to interact with the electrostatic surface potential of Phe481. (a-e) GVs of Phe and fluorinated 
Phe derivatives incorporated at position 481 on the background of Lys380Cys (a), Arg387Gln (b), Lys390Gln (c), Arg394Gln (d) and Arg377Lys (e) 
(n = 3-11 each, see also Supplementary Table S1); (f) AZFV 1/2 values plotted against the number of fluorines for Phe and fluorinated Phe derivatives 
incorporated at position 481 on the background of Lys380Cys, Arg387Gln, Lys390Gln, Arg394Gln and Arg377Lys. Lines indicate linear fits obtained 
with f=y0 + a*x, with y0 = 0 (mutants in black, WT in red). All data = mean ± s.e.m. 



initio calculations to further characterize such a possible repulsive 
electrostatic interaction. First, we aligned a negatively charged 
acetate ion with the face of a Phe and calculated the resulting 
energy over a range of distances. Not surprisingly, the energy 
required to maintain a given distance rapidly dissipates with 
increasing separation (Supplementary Fig. Sll). Plotting the 
experimentally determined free energy differences based on the 
fluorination of Phe48 1 against the predicted energies at an acet- 
ate-Phe distance of 7 A (similar to that observed in the Kvl.2/ 
Kv2.1 chimera open channel crystal structure 3 ), we observed a 
very strong inverse correlation (R = 0.98), with a slope of ( — )1.0 
(Supplementary Fig. 11). However, it should be noted that the 
(unknown) dielectric constant in the space between Phe481 and 
Glu395 would modulate the strength of the interaction, 
complicating a direct conclusion about the precise distance. 
Finally, we sought to test two possible relative orientations of the 
aromatic ring and the carboxyl group by calculating the 
interaction energies between a carboxyl group and Phe or 



fluorinated Phe derivatives in both an en face and an en edge 
orientation (Supplementary Fig. 11). Although the orientation of 
Phe481 and Glu395 varies among different structural and 
computational models 3,42,43 , the results from the ab initio 
calculations suggest that the experimentally observed data are 
consistent with an en face orientation, as the data for the en edge 
orientation predicts large positive and non-linear energy changes 
with fluorination, in contrast to the linear and negative changes 
observed experimentally. It should be noted that in principle, 
Phe481 could adopt different ro tamer conformations, but while 
we cannot rule out this possibility, it seems unlikely given the 
limited space around Phe481 (due to the close proximity of 
Tyr485, Phe402, Thr388 and Leu472). 

Discussion 

Our results show that Phe481, located at the distal end of S6, is 
unique among the numerous aromatic side chains found in this 
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area, in that fluorination induces a highly specific stabilization of 
the open state. The phenotype elicited through fluorination 
suggests that the negative electrostatic surface potential on the face 
of Phe48 1 is crucial for the observed effect. Furthermore, we used 
SCA to identify potential charged interaction partners near 
Phe481. Although four basic side chains were found to be close 
to Phe481, all displayed low SCA values and none altered the 
effect of fluorination at position 48 1 when neutralized. In contrast, 
neutralization of Glu395, a highly co-evolved acidic side in very 
close proximity to Phe481 (although located on the adjacent 
subunit), abolished the effects of fluorination at Phe481, raising 
the possibility of a repulsive interaction between the face of 
Phe481 and the carboxyl moiety of Glu395 in the channel open 
state (Fig. 7). This notion is not only consistent with the unusually 
high coupling score for the Phe481-Glu395 pair but further 
supported by ab initio calculations. It should be noted, however, 
that our data do not fully exclude the possibility of an indirect 
interaction between Phe481 and Glu395, such as would arise from 
another side chain located between Phe481 and Glu395 (with 
Tyr485 being a potential candidate, see Supplementary Fig. S10), a 
lipid head group located between the two side chains (an acyl 
chain is located within about 8 A of Phe481 in the Kvl.2/Kv2.1 
chimera open channel crystal structure, although no head group 
was resolved 3 ) or an indirect steric effect/local rearrangement 
caused by mutations at position 395. 

Our finding that Phe481 is highly sensitive to replacement with 
both naturally occurring and unnatural side chains is entirely 
consistent with previous studies on Shaker 12,44 and other Kv 
channels 15 ' 20 , which have shown that conventional mutations 
have potent functional effects and often result in nonfunctional 
channels. Further, it is noteworthy that the Phe481 equivalent 
aromatic side chain is absolutely conserved among all sodium 



Closed 




Figure 7 | Cartoon of a putative repulsive electrostatic interaction 
between Phe481 and Glu395 in the channel open state. Cartoon of the 
pore forming helices of two opposing potassium channel subunits in the 
closed (upper panel) and open (lower panel) conformation, respectively. 
VSD stands for voltage-sensor domain. The close proximity of the Glu395 
carboxyl and the Phe481 aromatic moiety in the open state are highlighted 
by the red oval; the electrostatic repulsion between the two is symbolized 
by the red arrows. The blue dots represent potassium ions along the 
permeation pathway. Note that for reasons of simplicity the model depicts 
an intrasubunit interaction between Glu395 and Phe481, while in the full- 
length channel the shortest distance between Glu395 and Phe481 is 
between adjacent subunits, thus indicating an intersubunit interaction. 
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channel isoforms 45 . Interestingly, an erythromelalgia- associated 
Val substitution of the equivalent Phe in domain III of the 
sodium channel isoform 1.7 (Navl.7) also results in a pronounced 
left-shift in the GV (and virtually no effect on activation 
kinetics) 45 ' 46 . These data may suggest a stabilization of the 
channel open state similar to that observed with fluorinated Phe 
derivatives at position 481 in Shaker. Reminiscent of studies 
highlighting the critical importance of Phe481, Glu395 has also 
been implicated previously as a determinant of pore stability: 
substitutions with the small, hydrophobic side chain Ala 
drastically left- shift the GV, similar to effects of trifluorination 
at Phe481 (Supplementary Table SI) 7 ' 19 , while replacement with 
Trp causes an extreme (>200mV) right-shift in the GV 11 , 
suggesting that charge is not the only critical parameter for 
channel function at position 395 (consistent with the observation 
that the Glu395Asp mutation causes drastic functional conse- 
quences and is a possible cause of episodic ataxia 47 ). Finally, 
structural studies have demonstrated the close physical proximity 
of the equivalent side chains to Phe481 and Glu395 in Kv 
channels 3,42 and, more importantly, comparison of open and 
closed structures of the full-length prokaryotic KcsA potassium 
channel show that the equivalent side chains to Phe481 and 
Glu395 approach each other by about 2 A in the channel open 
state, coming as close as 4 A (refs 48,49), a notion further 
supported by a recent all- atom simulation of Kv channel gating 
transitions 9 . 

It should be noted that aromatic molecules possess the unique 
capability to engage in attractive interactions with both 
cations (cation-pi interactions 24 ' 25 ) and anions (anion-pi 
interactions 44,50 ). However, in the present study, both of these 
possibilities seem very unlikely: if Phe481 was involved in a 
cation-pi interaction, fluorination would abolish this interaction 
and thus would not manifest as a stabilizing a certain state 
(rather, fluorination would destabilize a particular state, a 
phenotype that is not observed experimentally). Anion-pi 
interactions are also attractive in nature and result from heavily 
substituted aromatic systems to generate a positive quadrupole 
moment sufficient in strength to engage in attractive interactions 
with anions. These conditions are clearly not fulfilled by F 3 -Phe 
and its quadrupole moment near zero 50,51 . 

Is it possible that the dipole elicited through fluorination is 
responsible for the drastic stabilization of the channel open state? 
We believe this is unlikely for four reasons. First, the F-Phe and 
F 2 -Phe side chains both have dipole moments of around 1.9 debye 
and should thus result in an identical functional effect on the 
channel, which is not seen experimentally. Second, if the altered 
dipole of a fluorinated Phe derivative is responsible for the open- 
state stabilization, then mutations of other residues would be 
unlikely to eliminate the effect of fluorination at position 481. 
However, neutralizing Glu395 effectively abolished the effect of 
fluorination at Phe481. Third, the dipole moment of the Tyr side 
chain is comparable to that of F-Phe and F 2 -Phe (1.7 debye versus 
1.9 debye) but the functional effects vary drastically: although Tyr 
in position 481 leads to a dramatic (100 mV) right-shift in the 
GV, F-Phe and F 2 -Phe both left-shift the GV by 6 and 12 mV, 
respectively. Lastly, ab initio calculations, which factor in effects 
of fluorination on the aromatic dipole, show that an en edge 
interaction predicts an energetically non-linear effect, 
which also was not borne out in the functional consequence 
of fluorination at Phe481. Similarly, it is improbable 
that the observed effect of open-state stabilization is caused by 
a steric effect of the fluorinated Phe derivatives because the 
addition of a single fluorine adds 5 A 3 to the volume of a Phe side 
chain, while introducing a hydroxyl group adds 8 A 3 (www. 
molinspiration.com) but both have opposing functional 
effects (see above). 



Taken together, the data demonstrate that the electrostatic 
surface potential of a highly conserved Phe481 regulates open- 
state stability in Shaker channels and further support the notion 
of a repulsive interaction between Phe481 and Glu395, which 
could serve to inherently destabilize the channel open state in WT 
channels. Regardless of whether the interaction between Phe481 
and Glu395 is direct or indirect (modulated possibly through the 
Tyr485 hydroxyl), we speculate that upon channel opening, the 
negatively charged carboxyl moiety of Glu395 and the electro- 
negative surface potential on the face of Phe481 come in closer 
proximity than in the closed state and thus lead to an electrostatic 
repulsion that modulates open-state stability. Although rare, 
glutamate residues buried in hydrophobic pockets have been 
reported before 52-54 and the resulting energetic penalty has been 
suggested to be important for protein function 53 ' . The proposed 
mechanism of state-dependent electrostatic repulsion could be 
functionally relevant in other ion channel types, and may be 
a broadly employed energetic contribution in balancing 
conformational dynamics in other protein types as part of the 
general tradeoff between stability and function 56 . 

Methods 

Molecular biology and in vivo nonsense suppression. Residues 6-46 were 
removed from Shaker cDNA (in pBSTA) to abolish N-type inactivation. All con- 
structs contained the point mutations C301S, C308S and T449V (unless stated 
otherwise) and successful incorporation of these and all additional mutations was 
confirmed by automated sequencing. cRNA was transcribed from template cDNA 
with the mMessage mMachine kit (Ambion, Austin, TX, USA). Stage V-VI 
Xenopus laevis oocytes were prepared, and injected as previously described 57 . 
Oocytes were incubated for 8-48 h after injection (at 18 °C). Fluorinated Phe 
derivatives (4-F-Phe (F-Phe); 3,5-F 2 -Phe (F 2 -Phe) and 3,4,5-F 3 -Phe (F 3 -Phe)) were 
purchased from Asis Chem (Watertown, MA) and Sigma Aldrich (St Louis, MO). 
The in vivo nonsense suppression methodology has been described in detail 
elsewhere 22 . In brief, unnatural amino acids (aa) were protected with 
nitroveratryloxycarbonyl and were subsequently activated as the cyanomethyl ester. 
The resulting product was then coupled to the dinucleotide dCA (Dharmacon, 
Lafayette, CO) and ligated to a modified (G73) Tetrahymena thermophila tRNA. 
Immediately before co- injection with the channel cRNA, UV irradiation was used 
to deprotect the aminoacylated tRNA-aa from the nitroveratryloxycarbonyl 
protection group. Typically, 10-80 ng of tRNA-aa and 25-50 ng of channel cRNA 
were injected in a 50 nl volume. To test for nonspecific incorporation of 
endogenous amino acids and/or recharging of the injected tRNA, the cRNA 
together with a tRNA coupled to pdCpA (without an appended aa) were injected. 
Even at large depolarizations ( + 50mV) currents for Phe401TAG, Phe402AG, 
Phe481TAG, Phe484TAG and Tyr485TAG constructs were never larger than those 
recorded from oocytes injected with water, (<500nA, n = 4 each). 

Electrophysiology and fluorometry. Two electrode voltage-clamp was used to 
voltage- clamp potassium currents in standard Ringers solution (in mM): 116 NaCl, 
2 KC1, 1 MgCl 2 , 0.5 CaCl 2 , 5 HEPES (pH 7.4) using an OC-725C voltage clamp 
(Warner, Hamden, CT). Glass microelectrodes were backfilled with 3M KC1 and 
had resistances of 0.1-1 Mfl Typically, the holding potential was — 80 mV and 
depolarizing voltage-pulses were given for 50-500 ms in 5 and 20 mV increments 
for conventional two electrode voltage-clamp and VCF, respectively. For VCF 
experiments, oocytes were labelled on ice in a depolarizing solution containing 
(in mM) 100 KC1, 1.5 MgCl 2 , 0.5 CaCl 2 , 10 HEPES (pH 7.4) and 10 uM tetra- 
methylrhodamine-5-maleimide (tetramethylrhodamine-5-maleimide; Invitrogen, 
Carlsbad, CA) for 20-30 min. Oocytes were subsequently washed and stored in 
Ringers solution for up to 2 h on ice. The VCF set up consisted of an Olympus 1X5 1 
inverted microscope, a PMT400 photomultiplier tube (IonOptix, Milton, MA), a X- 
Cite 120 Q light source (Lumen Dynamics Group, Mississauga, ON, Canada) and a 
Lambda SC shutter to minimize photo bleaching (Sutter Instruments, Novato, CA). 
The fluorescence signal was acquired with a sampling frequency of 1 kHz and 20- 
50 sweeps were averaged for analysis; for display, fluorescence traces were offline 
filtered at 300-800 Hz. Conductance-voltage (GV) relationships were obtained by 
plotting the isochronal tail current amplitudes against the depolarizing pulse 
potential or by dividing the measured current amplitude at a given voltage by the 
driving force (assuming a potassium reversal potential of — 84 mV). The peak 
fluorescence signals were plotted as a function of the depolarizing pulse potential to 
obtain fluorescence-voltage relationships. All data = mean ± s.e.m. 

Kinetic modelling. The Zagotta-Hoshi-Aldrich model of Shaker gating was 
employed to model the effects of manipulating different rate constants 9 . 
IonChannelLab 1.0.5.3 was used to run the model and to simulate both steady-state 
curves and deactivation kinetics 58 . Deactivation time constants of the simulated 
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data were fitted with single exponentials in Clampfit 10.2.0.12 (Molecular Devices, 
Sunnyvale, CA). 

SCA and ab initio computations. SCA 35 was performed based on a multiple 
sequence alignment (MSA) representing 360 members of the Kv family, as used 
previously 37 . We used multiple iterations of a PSI-BLAST (e-score < 0.001) on a 
non-redundant database to generate the MSA, and used KvAP, Shaker and BK 
channel sequences for the initial query 59 . Only the sequence regions of the voltage- 
sensor and/or pore regions were used; further details of the MSA are described 
elsewhere 37 . SCA was performed using the SCA 5 MATLAB toolbox downloaded 
from the Ranganathan Lab website (http://systems.swmed.edu/rr_lab), deployed in 
MATLAB Ver. 7.10 (Math Works, Natick, MA). Instead of truncating the MSA 
based on gap -frequency, we used the whole alignment and then refined the data to 
only include those positions found on the Shaker sequence for analysis and figures. 
Ab initio computations gas phase Hartree-Fock 6-3 1G calculations were performed 
using Gaussian09/GaussView 5.0 (Gaussian, Wallingford, CT) as described 
elsewhere 57 . 
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